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A method is presented which corrects the one-fluid conformal solution viscosity model for size and mass 
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1. Introduction 

In a series of recent reports [1,2, 3], 1 Ely and Hanley have 
proposed a corresponding states method for the prediction 
of the viscosity and thermal conductivity of pure hydrocar- 
bons and their mixtures. This work was an extension of the 
previous work of Hanley [4, 5], which dealt with the trans- 
port properties of liquefied natural gas mixtures, to molecu- 
lar weight ranges corresponding to C 20 and other chemical 
types (e.g., aromatics). The method is based on a one-fluid, 
conformal solution concept and requires only pure compo- 
nent, equilibrium parameters such as the critical parame- 
ters as input. No transport data are required. 

Extensive comparisons of the predictions of the model 
with experimental data have been reported and are sum- 
marized in tables 1 and 2. In general the results are ex- 
cellent with the average absolute error between experiment 
and prediction being less than 8 percent for both pure 
fluids and mixtures. It was noted, however, that when the 
size difference of two binary mixture species becomes large 
(e.g., VyV c 2 ~ 6), the predictions of the viscosity model 
become markedly worse. This failure of the one-fluid model 
for viscosity has been explained by the nonequilibrium 
molecular dynamics studies of Hanley and Evans [6, 7]. 
These studies have shown that for mixtures of molecules of 
substantially different size, the mean density approximation 
inherent in the one-fluid theory for the binary pair distribu- 
tion function fails, even for a conformal mixture. A conse- 
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quence of this failure is that the local or ambient concentra- 
tion of the mixture components is not the same as the bulk 
concentration. This is shown in figure 1 for a 50/50 mixture 
of a conformal system whose size difference is two [7]. As 
one can see, the concentration of the larger component 
about a central large molecule (x 22 ) is greater than the bulk 



TABLE 1. Summary of One-Fluid Corresponding States Viscosity Results 
for Pure Fluids' 



Fluid Type 


N 


AAD 


BIAS 


n-Paraffins 


1301 


4.89 


-0.48 


i-Paraffins 


155 


21.17 


-21.17 


Alkenes 


58 


11.29 


7.85 


Cycloalkanes 


89 


40.56 


-40.56 


Alkylbenzenes 


155 


8.45 


-0.69 


Carbon Dioxide 


111 


4.75 


-4.53 


Overall 


1869 


8.42 


-4.10 



* AAD = Average absolute percent deviation. BIAS = Average percent 
deviation. 

TABLE 2. Summary of One-Fluid Corresponding States Viscosity Results 
for Binary Mixtures.' 



Mixture Type 


N 


AAD 


BIAS 


Alkane/Alkane 


303 


5.89 


-1.79 


Alkane/Cycloalkane 


24 


17.31 


-16.51 


Alkane/Alkylbenzene 


128 


7.41 


-0.01 


Overall 


455 


6.95 


-2.07 



* AAD = Average absolute percent deviation. BIAS = Average percent 
deviation. 

These results were obtained using the empirical size difference correc- 
tion proposed by Ely and Hanley [1]. 
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FlCURE 1. Local mole fractions for a conformal mixture of soft spheres 
with a mass ratio of 10 and size ratio of 2 plotted versus reduced inter- 
molecular separation r* [7]. 

concentration until one exceeds four or five molecular 
diameters. Since the viscosity (as well as other thermophy- 
sical properties) are determined from relatively short-range 
forces, the large component dominates the value of the mix- 
ture viscosity. Ely and Hanley [1] attempted to correct for 
this effect with an empirical relation based on the size ratios 
in the mixture. Although this function was somewhat suc- 
cessful (as is reflected in table 2), the size difference effects 
persist in the model predictions. 

In this manuscript, a systematic correction to the one- 
fluid viscosity model is proposed for size and mass dif- 
ference effects. This correction is based on the exact solu- 
tion of the Enskog model for a multicomponent mixture of 
hard spheres [8]. This approach has a rigorous foundation 
in the perturbation expansion of an equilibrium property of 
a fluid [9], but is empirical as applied to transport 
phenomena. In spite of this empiricism, the proposed cor- 
rection does improve the viscosity predictions for mixtures 
which exhibit large size and mass differences, for both the 
dense liquid and dilute gas states. 

Section 2 of this article summarizes the assumptions and 
working equations of the one-fluid, conformal solution 
viscosity model (CSVM). Section 3 discusses the hard sphere 
expansion model and describes the Enskog solution for a 
multicomponent mixture of hard spheres which is the ana- 



lytical formulation used to correct the CSVM. In section 4 
the predictions of the corrected and uncorrected models are 
compared with experimental data for both the dilute gas 
and high density fluids. Unfortunately, for methane/n- 
decane like systems where the size and mass difference ef- 
fects would be the most pronounced, no dilute gas experi- 
mental viscosities have been measured. For this reason, the 
model predictions are also compared to calculated Lennard- 
Jones viscosities. 

2. One-Fluid Viscosity Model 

In the one-fluid conformal solution viscosity model there 
are three basic assumptions: (1) the viscosity (rj) of a mixture 
at a density q, temperature T and composition [x a ] can be 
equated to the viscosity of a hypothetical pure fluid, i.e., 
Vmix (q, T, [x a )) =rf x (g, T); (2) the viscosity of the hypotheti- 
cal pure fluid may be evaluated via a corresponding states 
principle 



Vx(q,T) = t] (q ,T ) F„ 



(1) 



where i% is a dimensional factor defined below and (3) the 
reference fluid density and temperature (q and T„) may be 
evaluated via an extended equilibrium corresponding states 
principle [10] viz. 

Qo = Qh x ™dT =T/f x (2) 

where h x and f x are defined by the relations 

A^Q,T)=f x A^ Q h x , T/f x ) 

and 

Z x ( e ,T)=Z ( Q h x , T/f x ). 

A R denotes the residual Helmholtz free energy of the 
hypothetical or reference fluid (subscripts x or o, respective- 
ly) and Z is the compressibility factor, p/gRT. In eq (1), F, is 
given by 



fA-§-)'" f." v 



(3) 



where M denotes the mass. 

In order to apply the model to pure fluids or mixtures, 
analytical expressions for f„ h„ and M x as well as for the 
reference fluid equation of state and viscosity surface are 
required. In our previous work, methane was chosen as the 
reference fluid owing to the availability of pVT and viscosi- 
ty data for that fluid. The appropriate correlations have 
been reported previously [1] and will not be repeated here. 
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For a mixture, f„ h x , and M x must be obtained via mixing 
and combining rules for the corresponding mixture compo- 
nent parameters. We have adopted a set of one-fluid mixing 
rules given by the following 



f x = h: 1 £ E xjc^h^ 



h x = L E XaX h a0 



and 



M x = [L ExjvnlFfU 2 hV/Yr^K 



(4) 



(5) 



(6) 



The combining rules for the binary pair parameters (as 
denoted by an "a/3" subscript) are given by 



and 



U =(W 2 (i-U 
h a0 =\{hl"+hl"Y{\-U 



m a0 = 2m a m l[mi3 + m g ]. 



(7) 
(8) 

(9) 



In eqs (7) and (8) k a0 and l a0 are the binary interaction 
parameters which can be set equal to zero in viscosity pre- 
dictions. The parameters f Q and h g are the equivalent 
substance reducing ratios for the energy and volume for 
component a in the mixture. They are given by 



/- = {Tim) 6(v:, t;, Wo ) 



(10) 

(ii) 



where the subscript "c" indicates a critical value, "*" 
denotes a value reduced by the critical point and w is 
Pitzer's acentric factor. 6 and (f) are the shape factors of 
Leach and Leland [11, 12] whose detailed functional forms 
are given in reference [1]. T denotes the absolute tempera- 
ture, V is the molar volume and m is the mass. 

The mass mixing rule given by eq (6) was derived by 
Evans and Hanley [6] in their study of the viscosity of a mix- 
ture of conformal soft spheres. It arises by examining the 
potential contribution to the pressure tensor in terms of the 
nonequilibrium radial distribution function and thus is a 
mixing rule for the potential or in practice, high density 
contribution to the viscosity. This rule was adopted for all 
densities, however, since the emphasis of our previous work 
was on the dense fluid states. One might expect, therefore, 
that the CSVM might be somewhat less accurate for the 
dilute gas, kinetic regime where the mixture mass depend- 
ence is effectively proportional to m 1 ' 2 / 1 ' 2 h~ 2 ' 3 rather than 



m u2 f n h 4 ' 3 as given in eq (6). This point will be discussed 
further in section 4. 



3. The Enskog Correction 

Mansoori and Leland and their co-workers [9, 13] have 
proposed a conformal solution model for equilibrium ther- 
modynamic properties in which a dimensionless or reduced 
property of a mixture is expanded about the corresponding 
property of a hard sphere mixture. For example, if X mix 
represents the value of the real mixture property, one 
obtains 

X mix (Q,T,ix a }) = X»U{Qoi}, ix a })-XH s (eoi) + X {Q^T^ (12) 

In this equation X™ x denotes the value of the property in a 
mixture of hard spheres of diameters [a a ], [x a ] denotes the 
mixture composition, X" s is the property value in a pure 
hard sphere fluid of effective diameter a x (e.g., a one-fluid 
approximation) and X is the value obtained from a real, 
pure fluid reference substance, evaluated at the state point 
(q<»T ) where q^qo^Io 3 , and T' = TeJe x . In terms of inter- 
molecular potentials the parameters o and e correspond to 
the points where u(o)=0 and min(u)= — e. Both a x and e x 
are one-fluid parameters which must be evaluated via mix- 
ing rules. The difficulty in applying this approach lies in 
assigning values to a and e for the mixture components. 
One possible approach is to assume that o 3 ~ V c and e ~ 
T c which leads to factors such as those given in eqs (10-1 1). 
The choice of parameters used in this work will be discussed 
in section 4. 

Although transport properties cannot, in general, be ex- 
panded in a perturbation series [14], it is tempting to apply 
the hard sphere expansion (HSE) formalism to the confor- 
mal solution viscosity model presented in section 2. For- 
mally this may be written as 



VUQX M, {m a }) = v^iiQO 3 ,}, M, {m a }) - 



tf'iQO 3 ,, m x ) + Ueo, T ) F, 



(13) 



where the notation is the same as defined previously. Note 
that in the case of transport we must also consider the 
masses of the particles {m a }. 

In practice, we do not have an exact model for the viscos- 
ity of a hard sphere fluid (pure or mixed) at all densities. For 
this reason, the Enskog model [15], which has been solved 
for a multicomponent mixture of hard spheres by Tham and 
Gubbins [8], was selected to calculate rf" x and r)" 5 . Their 
solution is given by 

r?™ = E ft(r, G ) Y,+^LL x iXj u°j y u 
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where (12-6) intermolecular potential value given by Reid, et al. 

[16]. This choice, although resonable, is still somewhat arbi- 
trary. The mixing rules for the one-fluid hard sphere system 
were chosen to be consistent with those used in the CVSM, 

and the 0, are the solutions to the set of linear equations e 9 s ( 5 > and ( 6 )' V1Z ' 

defined by the following _ _ , 

of = E L, x a x b a% b 



Y,=x,[l +%XxjM Jiy(/ ] 



(16) 



ZB u Pj = Y, 



where 



B, = 2 E x,xmt Mi (*JT 1 [d + % ^) «<, - %(%) 8 J( ] . 

In these equations, fe y = 27ro-£,/3, a y = x Mp t + o>), Mj, = 
mjl{m, + mj), u%- b^nifo, rfi, = 5 / 16 (^4 m y = 
2mtMj,, y u = n byg^aj), n = number density, g y is the 
hard sphere radial distribution function for the ij pair, and 
6 is the Kronecker delta function. Although these equations 
are somewhat complex, they may be readily solved on a digi- 
tal computer. 

The final model, which we shall call the hard sphere 
expansion-conformal solution viscosity model (HSE-CSVM), 
is given by 

*-«. {Q,T, M, {m a }) = At? ««» + ,. ( Do , T.) F n (14) 
with 

„ENSKOG _ -ENSKOG . 



±r, ENSKOC = V E m T OC ({Q<Z}. M, M) ~ 



(qoJ, m x ) 



(15) 



4. Results 

In order to apply eqs (14) and (15) values for the hard 
sphere diameters a a must be chosen and mixing rules for 
the one-fluid values a x and m x in the hard sphere system 
must be selected. As was mentioned previously, consid- 
erable freedom exists for the choice of the molecular diam- 
eters. Unfortunately, the hard sphere contributions to the 
viscosity are rather sensitive to this choice, thus several 
possibilities were considered. The relationship which was 
selected is given by 



o. = (^A o /3.058W o )w 



(18) 



where h„ is the ratio defined in eq (1 1) and is obtained in the 
CSVM calculations and /V„ is Avogadro's number. The fac- 
tor 3.058 was chosen so that methane would have a diameter 
of 3.758 10" 10 m which corresponds to the Lennard-Jones 



and 



M x =[L12x a x b m^ai b ] 2 a- x 



where a ab = Vz {a a + o b ) and m ab is defined in eq (9). 



4.1 Results for the Dilute Gas Limit 

The correction was first tested on the dilute gas viscosity 
of a methane/propane mixture with the results being given 
in table 3. This table gives the experimental data [17], 
CSVM predictions and the HSE-CSVM predictions ob- 
tained using eqs (14-16). As one can see, the HSE-CSVM is 
consistently more accurate with average absolute percent- 
age deviation being 3.6 percent as compared to 6.29 percent 
for the CSVM. 

The size and mass difference in the methane/propane 
system is not very great— a, lo 2 ~ 1.4 and m l lm 2 — 3. One 
would expect that the effect of the correction would be more 
pronounced in a system like methane/n-decane where the 
size and mass ratios are 1.8 and 9, respectively. Unfor- 
tunately no experimental measurements for the dilute gas 
viscosity of this mixture have been reported. For this 
reason, the Lennard-Jones gas viscosities of this system 
were calculated using the standard kinetic theory formalism 
[15]. Although it is impossible to assess the absolute accu- 
racy of these values, they do serve as a rational basis upon 
which the HSE-CSVM and CSVM may be compared. 

In order to perform the dilute gas Lennard-Jones (LJ) 
calculations, values for the intermolecular potential 
minimum e a for the mixture components must be chosen as 
well as values for the ct„. The o a were obtained from eq (18) 
and ejk where k is Boltzmann's constant were calculated 
from the empirical relation 

ejk = n/Jl.282 

where/, is the CSVM reducing parameter given in eq (10). 
The factor 1.282 was chosen so that methane would have a 
value of dk of 148.6 which is the Lennard-Jones value 
reported by Reid, et al. [16]. 

Table 4 summarizes the results obtained with the LJ 
kinetic theory model and the HSE-CSVM and CSVM at 
three compositions and four temperatures. Comparing the 
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TABLE 3. Comparison of Calculated and Experimental [17] Dilute Gas Viscosities of Methane/Propane Mixtures at 1 atm. 



Composition 


Temperature 
K 


7j„ p , 10 7 Pa.s - 




-CSVM 
Vcolc 


HSE-CSVM 




mol %, methane 


10 7 Pa«s 


% 


W Pa.s 


% 


22.07 


310.928 


89.2 


93.1 


4.37 


90.6 


1.57 




344.261 


98.8 


103.2 


4.45 


100.7 


1.92 




377.594 


108.3 


113.5 


4.80 


111.0 


2.49 




410.928 


118.0 


123.8 


4.92 


121.3 


2.80 


38.78 


310.928 


93.0 


99.2 


6.67 


95.6 


2.80 




344.261 


103.0 


109.9 


6.70 


106.5 


3.40 




377.594 


113.0 


120.6 


6.73 


117.0 


3.54 




410.928 


122.5 


131.2 


7.10 


127.7 


4.24 


61.39 


310.928 


99.6 


107.3 


7.73 


103.6 


4.02 




344.261 


109.6 


118.4 


8.03 


114.9 


4.84 




377.594 


119.5 


129.4 


8.28 


125.9 


5.36 




410.928 


129.2 


140.1 


8.44 


136.7 


5.80 


79.10 


310.928 


107.2 


112.6 


5.04 


110.2 


2.80 




344.261 


117.4 


123.7 


5.37 


121.4 


3.41 




377.594 


126.8 


134.5 


6.07 


132.3 


4.34 




410.928 


136.6 


144.8 


6.00 


142.8 


4.54 




Average absolute percent deviation 




6.29 




3.62 


TABLE 4. Comparison 


of Calculated and Lennard-Jones (12-6) Dilute Gas Viscosities < 


if Methane/n-Decane 


Mixtures and CSVM Predictions 


Composition 


Temperature 
K 


t,£,10 7 Pa-s - 




„CSVM 


-HSE-CSVM 




mol %, methane 


10 7 Pa«s 


% 


10 7 Pa.s 


% 


25.0 


300 


46.5 


54.9 


18.06 


43.4 


-6.67 




400 


65.4 


76.2 


16.51 


62.8 


-3.98 




500 


86.4 


99.2 


14.81 


84.1 


-2.66 




600 


108.8 


124.2 


14.15 


107.6 


-1.10 


50.0 


300 


56.4 


74.9 


32.80 


50.1 


-11.17 




400 


79.2 


103.1 


30.18 


74.9 


-5.43 




500 


103.7 


133.4 


28.64 


102.4 


-1.25 




600 


129.1 


165.2 


27.96 


132.3 


2.48 


75.0 


300 


75.2 


104.2 


38.56 


70.6 


-6.12 




400 


103.6 


140.5 


35.62 


104.8 


1.16 




500 


132.2 


177.0 


33.89 


140.5 


6.28 




600 


160.7 


211.5 


31.61 


175.4 


9.15 




Average absolute 


1 percent deviation 




20.17 




3.59 




Average percent 


deviation 




20.17 




-2.40 



results we see a substantial difference in the CSVM predic- 
tions and the LJ calculations. In the case where the methane 
content is high, (75 percent) the HSE-CSVM is in better 
agreement with the LJ calculations by 20-30 percent. Also, 
the overall absolute percentage deviation for the HSE- 
CSVM is 3.6 percent as opposed to 20 percent for the 
CSVM version. Although the accuracy of the LJ calculations 
is uncertain, similar calculations on the methane/propane 
system where experimental data do exist yielded agreement 
between the LJ and experiment of better than 2 percent. 
One can postulate, therefore, that discrepancies similar to 
those shown in table 4 would exist between the CSVM and 
real experimental data for the methane/n-decane system. 

Note that in both comparisons the CSVM predicts values 
which are too large in the dilute gas limit. We attribute this 



to the potential or high density mass mixing rule which is 
used in F, in eq (3). Thus, we see that in the low density 
limit the HSE Enskog correction is negative and lowers the 
predicted viscosity. 



4.2 High Density Results 

The initial motivation for this work was to obtain a cor- 
rection for the effect of size differences on the predicted 
high density viscosity of systems like methane/n-decane mix- 
tures. Figure 2 compares the experimental [18] and CSVM 
predicted viscosities of methane/n-decane mixtures as a 
function of reduced density at three different compositions. 
Note that the predictions are worst for the high methane 
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FicliRE 2. Comparison of calculated and experimental viscosities of 
methane/n-decane mixtures using the uncorrected model. Note that the 
predictions are worst (~ 30 percent error) for the highest methane concen- 
tration. 
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FIGURE 3. Comparison of calculated and experimental viscosities of 
methane/n-decane mixtures using the corrected model. Note that both the 
bias and average absolute deviations are improved. 



composition and improve with increasing decane content. 
This is somewhat surprising since methane is the reference 
fluid in the CVSM calculations. This result was explained in 
the introduction and is attributable to a failure of the one- 
fluid theory to adequately represent the size difference ef- 
fects in the high density region. 

Figure 3 compares the experimental and calculated 
results for the HSE-CSVM. We see in this case that there is 
a marked improvement in both the bias and overall devia- 
tion. Note that in this case the HSE correction is positive. 
The density dependence of the HSE correction is illustrated 
in figure 4 which shows -ty- = (r) E m N ™ OG - r)* NSKOG )IT* 
plotted against reduced density at a size ratio of o x la 2 = 2.0 
and mass ratio of m 1 lm 2 = 8 at three compositions. This 
corresponds approximately to a methane/decane like 
system. In figure 4 the density was reduced by an approx- 
imation to the critical density of the mixture given by Qc 1 = 
3.058 N (x,of + x 2 cr|). This figure demonstrates that the 
correction decreases in magnitude with increasing concen- 
tration of the larger component (x x ) and is small and 
negative below reduced densities of 1.5. Also the correction 
increases sharply above reduced densities of 1.5. 

Table 5 compares the overall predictions for both the cor- 
rected and uncorrected model for 24 binary systems. The 
literature references for the experimental data are given in 
[1]. The systems where the size difference is large [a 2 la 1 = 
{V c 2 /V\f>\ show substantial improvement while those of 
similar size are essentially unchanged. It may be possible to 
further improve the corrected model by a judicious choice 
of the o*„. This possibility is currently being examined. 
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TABLE 5. Summary of Calculated and Experimental Dense Fluid Binary Mixture Viscosities: 



Component 1 


Component 2 


aja t 


N 


AAD™ 


BIAS"™ 


A KT)HSE-CSVM 


BIAS""-"™ 


Methane 


Propane 


1.273 


134 


5.91 


-5.29 


4.62 


-3.45 




n-Nonane 


1.778 


32 


6.37 


-5.58 


4.12 


-2.61 




n-Decane 


1.839 


71 


14.43 


-14.43 


5.35 


-1.54 


2,3-Dimethylbutane 


n-Hexane 


1.014 


2 


5.32 


-5.32 


5.31 


-5.31 




n-Octane 


1.110 


2 


6.03 


-6.03 


5.65 


-5.65 


n-Hexane 


n-Tetradecane 


1.304 


10 


2.15 


-1.15 


1.92 


0.54 




n-Hexadecane 


1.356 


26 


4.04 


-3.85 


2.59 


-1.97 


n-Heptane 


n-Dodecane 


1.185 


3 


2.47 


2.47 


3.44 


3.44 




n-Tetradecane 


1.242 


3 


1.19 


0.10 


1.82 


1.51 




n-Hexadecane 


1.291 


3 


3.03 


-3.03 


2.52 


-1.32 




n-Octadecane 


1.333 


2 


2.95 


-2.95 


1.92 


-0.71 


n-Octane 


n-Decane 


1.074 


2 


3.09 


3.09 


3.27 


3.27 


n-Tetradecane 


n-Hexadecane 


1.040 


11 


2.31 


2.08 


2.32 


2.10 


Benzene 


n-Hexane 


1.130 


15 


5.85 


-2.70 


5.74 


-2.28 




n-Heptane 


1.186 


3 


4.68 


4.68 


5.85 


5.85 




2,2,4-Trimethylpentane 


1.218 


26 


13.14 


-13.14 


12.46 


-12.46 




n-Decane 


1.329 


3 


4.45 


0.73 


5.47 


3.46 




n-Dodecane 


1.406 


3 


2.87 


2.87 


6.47 


6.47 




n-Tetradecane 


1.473 


3 


2.96 


-1.47 


3.63 


2.55 




n-Hexadecane 


1.531 


3 


3.75 


-2.80 


3.76 


1.46 




n-Octadecane 


1.581 


3 


2.99 


-2.99 


2.73 


1.97 


Toluene 


n-Heptane 


1.110 


21 


5.15 


5.15 


5.32 


5.32 




n-Octane 


1.157 


20 


9.03 


9.03 


9.50 


9.50 




2,2,4-Trimethylpentane 


1.140 


28 


6.61 


-4.64 


6.62 


-4.37 




Overall 




429 


7.45 


-5.23 


5.31 


-1.86 



The CSVM results were obtained without the empirical size difference correction give in [1]. 



5. Summary and Conclusions 



6. References 



We have shown that a relatively simple correction to the 
one-fluid conformal solution viscosity model may be ob- 
tained from the Enskog hard sphere theory. This function 
effectively corrects for errors in the mass mixing rules at low 
density and also for size difference effects at high density. 
Even though the effective correction is not substantial for 
systems of similar size and mass, it does systematically cor- 
rect the model predictions as is reflected in the BIAS as 
given in table 4. In addition, when the size and mass dif- 
ference is large, as in the methane/n-decane system, the im- 
provement in the model predictions is substantial. Further 
work is in progress to extend this approach to other struc- 
tural features such as branching and polarity which are not 
adequately handled by the one-fluid corresponding states 
model. 
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